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Abstract 

We derive inequalities between the area, the angular momentum and the charges 
for axisymmetric closed outermost stably marginally outer trapped surfaces, embed- 
ded in dynamical and, in general, non-axisymmetric spacetimes satisfying the Einstein- 
Maxwell-dilaton-matter equations. In proving the inequalities we assume that the dila- 
ton potential is nonnegative and that the matter energy-momentum tensor satisfies the 
dominant energy condition. 

1 Introduction 

The dynamical black holes are a serious challenge to the present day investigations in general 
relativity and alternative theories of gravity. The black hole dynamics is very difficult to 
be studied within the framework of the existing theoretical scheme and consequently our 
understanding of the dynamical black holes is not so deep as for the isolated stationary black 
holes. In this situation, the derivation of certain estimates and inequalities on the physical 
characteristics of dynamical black holes based mainly on "first principles" and independent 
of the specific features of the dynamical processes, is very important. Within the general 
theory of relativity, lower bounds for the area of dynamical horizons in terms of their angular 
momentum or/and charge were given in 0]H|7]|, generalizing the similar inequalities for the 
stationary black holes lf8l- lfT0l . These remarkable inequalities are based solely on general 
assumptions and they hold for any axisymmetric but otherwise highly dynamical horizon in 
general relativity. For a nice review on the subject we refer the reader to ifTTI . 

A natural problem is to find similar inequalities in alternative theories of gravity which 
generalize Einstein theory. An example of such a theory is the so-called Einstein-Maxwell- 
dilaton gravity which naturally arises in the context the generalized scalar-tensor theories of 
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gravity, the low energy string theory lfT2l[T3l , Kaluza- Klein theory [fT4ll . as well as in some 
theories with gradient spacetime torsion |fT5l . 

The field equations of Einstein-Maxwell-dilaton gravity with matter are presented below 
in eq. ©■ A characteristic feature of this theory is the coupling between the scalar field 
(dilaton) (p and the electromagnetic field F ab and this coupling is governed by a parameter 
y (called dilaton coupling parameter). The static and stationary isolated black holes in 4D 
Einstein-Maxwell-dilaton theory were extensively studied in various aspects during the last 
two decades. The classification of the isolated stationary and asymptotically flat black holes 
in Einstein-Maxwell-dilaton gravity was given in [16] for dilaton coupling parameter y satis- 
fying < y 2 < 3. The sector of stationary Einstein-Maxwell-dilaton black holes with dilaton 
coupling parameter beyond the critical value y 2 = 3 is extremely difficult to be analyzed ana- 
lytically. Most probably the black hole uniqueness is violated in this sector as the numerical 
investigations imply ifTTll . 

In the present paper we derive some inequalities between the area, the angular momen- 
tum and the charges for dynamical black holes in Einstein-Maxwell-dilaton gravity with 
a non-negative dilaton potential and with a matter energy-momentum tensor satisfying the 
dominant energy condition. 

2 Basic notions and setting the problem 

Let (M,g) be a 4-dimensional spacetime satisfying the Einstein-Max well-dilaton-matter 
equations 



where g a b is the spacetime metric and V a is its Levi-Civita connection, G a b = Rab — ^Rgab 
is the Einstein tensor. F ab is the Maxwell tensor and J a is the current. The dilaton field is 
denoted by cp, V(cp) is its potential and y is the dilaton coupling parameter governing the 
coupling strength of the dilaton to the electromagnetic field. The matter energy-momentum 
tensor is T ab . We assume that T a b satisfies the dominant energy condition. Concerning the 
dilaton potential, we assume that it is non-negative (V (cp) > 0). 

Further we consider a closed orientable 2-dimensional spacelike surface ® smoothly 
embedded in the spacetime fM. The induced metric on ® and its Levi-Civita connection 
are denoted by q a t> and D a . In order to describe the extrinsic geometry of *B we introduce 
the normal outgoing and ingoing null vectors l a and k a with the normalization condition 
g(l,k) — l a k a = — 1. The extrinsic geometry then is characterized by the expansion & l , the 
shear o l ab and the normal fundamental form £l l a associated with the outgoing null normal l a 



Rab - \Rgab = 2V fl (pV,cp - g ab V c (pV c <p + 2e- 2 ^ (F ac F b c - 8 -fF cd F cd } 

-IV^gab + SnTab, 

V[ a F bc ] = 0, 



(1) 
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and defined as follows 



® l = q ab Vak, 



dab = tfalftcld - \® l qab, 



-k c q d a V d l c 



(2) 
(3) 
(4) 



In what follows we require H to be a marginally outer trapped surface (i.e. & 1 = 0) and 
(B to be stable (or spacetime stably outermost in more formal language) ll22~Tl - |[20lU 41. The 
last condition means that there exists an outgoing vector V a = X\l a — X2k a with functions 
A,i > and %2 > such that dy® 1 > 0, with 8y being the deformation operator on 1> [fT9l - 
11211 . In simple words the deformation operator describes the infinitesimal variations of the 
geometrical objects on 1> under an infinitesimal deformation of 1> along the flow of the vector 
V a . 

As an additional technical assumption we require 1> to be invariant under the action of 
U(l) group with a Killing generator r\ a . We assume that the Killing vector r\ a is normalized 
to have orbits with a period 2n. Also we require that (B is axisymmetrically stableQ and 
fr|/ a = £^k a = and £r\Cl' a = £y\F a b = £ri9 = 0, where F is the projection of the Maxwell 
2-form on r B. 

From the axisymmetric stability condition one can derive the following important in- 
equality valid for every axisymmetric function a on 2? 



|£>a|2 + Itf s a 2 



dS> 



a 2 \CP \l + apVl« + G ab al a (ak b + $l b ) 



dS, 



(5) 



where | . \ q is the norm with respect to the induced metric q a b, dS is the surface element 
measure on (B, Rrg is the scalar curvature of r B, Q."^ — T\ a £l l a and (3 = (xk\/%i- 
At this stage we can use the field equations (OQ) which gives 



|Da| 2 + ^? s a 2 



dS> ^|a 2 |^ 11 | 2 + a(3|o / | 2 + a 2 |D(p| 2 + 2a 2 V((p) 



(6) 



+2a|3(/ fl V fl (p) 2 + aV 2w [E\ +b\ ] + 2a^e~ 2 ^ (iiF) a (iiF) a + SnT ab al a (ak b + p/*)} dS 

where Ej_ = i^UF and B± = ijji-kF. All terms on the right hand side of the above inequal- 
ity are non-negative. Indeed, for the last term we have $nT a i,(xl a ((xk b + $l b ) > since the 
energy-momentum tensor of matter satisfies the dominant energy condition. We also have 
2a(3e~ 2w (z'/i 7 ) fl (//i 7 ) a > since the electromagnetic field satisfies the null energy condition 
and aP > 0. 

Considering now the inequality for a = 1 and applying the Gauss-Bonnet theorem^ we 
find that the Euler characteristic of (B satisfies 



Euler(-B) > 0, 



(7) 



H.e. axisymmetric and stable with axisymmetric functions A-i and X2 

2 

4tc(1 — g) where g is the genus of (B 



2 According to the Gauss-Bonnet theorem, the Euler characteristic is given by Euler^) — A f%R%dS 
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which shows that the topology of 2? is that of a 2-dimensional sphere S 2 . 

Discarding the following non-negative terms aPla'l 2 , 2a 2 V((p), 2a(3(/ a V ( . ! (p) 2 , 
2a$e- 2 W(iiF) a (iiF) a and %%T ab al a (ak b + $l b ) we obtain 



f |Da| 2 + ^a 2 dS> J^a 2 {\nX+\ D( ?\q + e ~ 2y(? [ E i+ B l}} ds ■ 



Proceeding further we write the induced metric on 2? in the form 

dl 2 = e 2C - a dQ 2 + e a sin 2 Qd<$> 2 , (9) 

where C is a constant. The absence of conical singularities requires o|e=o = c|e=jr = C. It is 
easy to see that the area of 2? is given by A — Ane c . Regarding the 1-form Q, l a , we may use 
the Hodge decomposition 

Q l = *dG) + dq, (10) 

where * is the Hodge dual on 2>, and (0 and q are regular axisymmetric functions on 2?. Then 
we obtain 

Q. y] =i n *d& (11) 

since q is axisymmetric and i-^dc, = £nC, = 0. 

We can also introduce electromagnetic potentials <J» and *F on 2? defined bjH 

d& = B ± *r\, (12) 
d¥ = e - 2y(? E ± *T\. (13) 

It turns out useful to introduce another potential % instead of co which is defined by 

dx = 2Xda-2<t>d x i' + 2 x i>d<i>, (14) 

where X = q a b^ ar \ b is the norm of the Killing field x\ a . This step is necessary in order to 
bring the functional I* [X A ] defined below, in the same formal form as in the stationary case. 
The electric charge Q and the magnetic charge P associated with 2? are defined as follows 

2=7-/ e- 2 ^ E± dS, (15) 
471 J<b 

P = — [ B ± dS. (16) 
4k J<b 

We also define the angular momentum J associated with 2? 

J=^T ( aV5+^- I (^e- 2 ^E ± -^B ± )dS, (17) 

o7t J<B oil J<B 

3 We denote the Killing vector field r| and its naturally corresponding 1-form with the same letter. 
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where the first integral is the contribution of the gravitational field, while the second integral 
is the contribution due to the electromagnetic field [fT6ll . 

Using the definitions of the potentials x ¥, $ and % one can show that the electric charge, 
the magnetic charge and the angular momentum are given by 

y(7c)-y(o) g(w)-g(o) xfr)-x(Q) fl81 

W 2 2 8 

Since the potentials *P, 4> and % are defined up to a constant, without loss of generality 
we put = - Y(0) = 2, 4>(ti) = -4>(0) = P and x(n) = -jc(0) = 47. 

Going back to the inequality ®, choosing a = e c ~°/ 2 , and after some algebra we obtain 

2(C+1) > ^^|a + ^-|Da| 2 + ^2|Dx + 24>7J>^-2 t I'D4>| 2 

+ I e -2Y9| Z )4,|2 + I e 2 Y 9| D »p|2 + |£>(p| 2 1 rfS , (19) 
A^ A^ J 

where the norm | . | and the surface element cISq are with respect to the standard usual round 
metric on S 2 . Taking into account that A = 4ne c the above inequality is transformed to the 
following inequality for the area 

W 7 ^- 2 )/ 2 , (20) 

where the functional /[X A ], with X A = (X,%,4>, l I / ,(p), is defined by the right hand side of 
Gil), i.e. 

I [ xA } = ^- [ |o + ^-|Do| 2 + -4t|Dx + 24>D^-2^D4>| 2 
2ll j<b y 4 4a z 

+ _L e -2Y<P| jD ^>|2 + 1 2 Y q)| D ^|2 + \ D( pA dSo . (21) 
XX J 

In order to bring the action into a form more suitable for the further investigation we 
express Do by the norm of the Killing field T] (i.e. e a = X / sin 2 0) and introduce a new 
independent variable X = cosG. In this way we obtain 



I[X A 



2 1 (d% ^d*¥ ^ Tl d®\ 2 

, A I — 4- 24> 2 V F — 

4X 2 V dx J 4X 2 \dx dx dx ) 



e -2v? fd®\ 2 + e 2 ™ fd*¥\ 2 + /dcp 



X \dx J X \dx J \dx 



1 



1-x 2 



dx. (22) 



At this stage we introduce the strictly positive definite metricB 



= (wrti- = ^ 2+te+ ^f- 2W< " )2 + £^!!+£!!^ + v (23) 

4a z a 



4 It is worth mentioning that the continuous rotational 0(2) symmetry in the case of Einstein-Maxwell 
gravity degenerates here to the discrete symmetry ±4> < — > ±1* and cp < — > —(p. 
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on the 5 -dimensional Riemannian manifold 9\f\ = {(X, %, 4>, *P, cp) G M 5 ;X > 0}. In terms of 
this metric the functional I[X A ] is written in the form 

Let us summarize the results obtained so far in the following 

Lemma. Let 'B be a smooth, spacetime stably outermost axisymmetric marginally outer 
trapped surface in a spacetime satisfying the Einstein-Maw ell- dilaton-matter equations ©. 
If the matter energy -momentum tensor satisfies the dominant energy condition and the dila- 
ton potential is non-negative, then the area of H satisfies the inequality 

A > W 7 ^- 2 )/ 2 , (25) 

where the functional I[X A ] is given by ([2?]) with a metric Gab defined by (|2J|). 

In order to put a tight lower bound for the area we should solve the variational problem for 
the minimum of the functional I[X A ] with appropriate boundary conditions if the minimum 
exists at all. Since the first two terms as well as the last term in I[X A ] are in fact boundary 
terms, the minimum of I[X A ] is determined by the minimum of the reduced functional 

In order to perform the minimizing procedure we have to specify in which class of func- 
tions X A — (XjX,^, 1 !',^), the functional is varied. From a physical point of view the 
relevant class of functions is specified by the natural requirements (%, <J>, cp) e C°°[— 1, 1], 

o = In (-f^i) G C°°[-l, 1] with boundary conditions 4>(x = -1) = -<&(% = 1)=P, = 

-1) = -V(x = 1) = Q and %(x = -1) = -%{% = 1) = 47. 

Lemma. For dilaton coupling parameter satisfying < y 2 < 3, there exists a unique smooth 
minimizer of the functional I*[X A ] with the prescribed boundary conditions. 

Proof. Let us consider the "truncated" functional 

/ T 2 dX A rlY B 

(l-x 2 )G AB ^^dx (27) 

with boundary conditions X A (x\ ) , X A (X2) for — 1 < X\ < X2 < 1 . Introducing the new variable 
t = \ In ({^) the truncated action takes the form 

rh dY A rlY B 

U[X A }[t 2 A} = I GAB^^dt, (28) 

which is just the geodesic functional in the Riemannian space (l7\£, Gab)- Consequently the 
critical points of our functional are geodesies in !7\£. However, it was shown in |fT6l that for 
< y 2 < 3 the Riemannian space (fA£, Gab) is simply connected, geodesically complete and 
with negative sectional curvature. Therefore for fixed points X A (t\) and X A {t2) there exists 
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a unique minimizing geodesic connecting these points. Hence we conclude that the global 
minimizer of Z* [X A ] fa, h] exists and is unique for < y 2 < 3. Since (fA£, Gab) is geodesically 
complete the global minimizer of 4[X A ][f2,?i] can be extended to a global minimizer of 

Even in the cases when the global minimizer of exists, there is another serious 

problem in Einstein-Maxwell-dilaton gravity. In Einstein-Maxwell gravity the lower bound 
for the area is clear from physical considerations and there is a completely explicit solu- 
tion realizing it, namely the extremal Kerr-Newman solution. So the approach is to formally 
prove that the area of the extremal Kerr-Newman solution is indeed the lower bound. The sit- 
uation in Einstein-Maxwell-gravity is rather different. Contrary to the Einstein-Maxwell case 
where the Euler-Lagrange equations can be solved explicitly, in Einstein-Maxwell-dilaton 
gravity the corresponding Euler-Lagrange equations are not integrable for general dilaton 
coupling parameter y. So there is no hope to find explicitly the sharp lower bound for the 
area in Einstein-Maxwell-dilaton gravity with arbitrary y. That is why our approach here 
should be different in comparison with the Einstein-Maxwell gravity. 

3 Area-angular momentum-charge inequality for critical 
dilaton coupling parameter 

The main result in this section is the next theorem: 

Theorem. Let Hbe a smooth, spacetime stably outermost axisymmetric marginally outer 
trapped surface in a spacetime satisfying the Einstein-Maw ell- dilaton-matter equations (0) 
with a dilaton coupling parameter y 2 = 3. If the matter energy -momentum tensor satisfies 
the dominant energy condition and the dilaton potential is non-negative, then the area of H 
satisfies the inequality 



JZ>Sn^\Q 2 P 2 -J 2 \, (29) 

where Q, P and J are the electric charge, the magnetic charge and the angular momentum 
associated with ( B. The equality is saturated only for the extremal stationary near horizon 
geometry of the y 2 = 3 Einstein-Maxwell-dilaton gravity with V((p) = and T a b = 0. 

Proof. For the critical coupling (fA£, Gab) is a symmetric space with a negative sectional 
curvature [16J. In fact fA£ is SL(3,R)/0(3) symmetric space and therefore its metric can be 
written in the form 

dh 2 = -Tr(M~ 1 dMM 1 dM) , (30) 
8 

where the matrix M is symmetric, positive definite and M £ SL(3,R). After tedious calcula- 
tions it can be shown that 



M = e 2 ^ 



/ x + 44>V 2 ^ /I( P +X- l (x + 2^) 2 le-^Q + lX-^x + l^)*? X- 1 (x + 24> v I , ) 
2 e - 2 ^^ + 2X- l (x + 2^) x i i e-^ + A^X- 1 mfX 1 
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In terms of the matrix M the functional I[X A ] becomes 



dx. (31) 



The Euler-Lagrange equations are then equivalent to the following matrix equation 

(l -x L )M i — = 0. 



dx 



dx I 



(32) 



Hence we find 



2 . dM 



(\-x z )M- l — = 2A 1 

CI I 



(33) 



where A is a constant matrix. From detM = l it follows that TrA = 0. Integrating further we 
obtain 



1 +x 

M = Mq exp I In A 

1 1 -x 



(34) 



where Mq is a constant matrix with the same properties as M and satisfying A t Mq = MqA. 
Since Mq is positive definite it can be written in the form Mq = BB T for some matrix B with 
| det5| = 1 and this presentation is up to an orthogonal matrix O (i.e. B — > BO). This freedom 
can be used to diagonalize the matrix B T AB T ~ 1 . So we can take B T AB T ~ 1 = diag(a i , «2, #3) 
and we obtain 



M = B 



( (sr \ 

(i±ir 
v eery 



5' 



(35) 



The eigenvalues a ; can be found by comparing the singular behavior of the left hand and right 
hand side of (1331 ) at x — > ±1. Doing so we find, up to renumbering, that a\ = 0, = — 1 and 
as = 1. The matrix 5 can be found by imposing the boundary conditions which gives 



/ 



B 



P 2 Qe 



(9++<P- 



^\P 2 Q 2 -J 2 \ 



Je 



(<P++9- 



-4?(<p++<p-) 



V 2 v /|/ >2 e 2 -/ 2 i 



(-2J + PQ)eV3 



(36) 



where 



<p± = <p(x = ±l), 
o p = km j 



g(x = ±1). 



(37) 
(38) 
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Taking into account that | det5| = 1 we find 

e°» = 2^\P 2 Q 2 -J 2 \. (39) 

Now we are ready to evaluate the minimum of the functional I[X A ] . Substituting (1331 in (1311) 
we see that the last two terms cancel each other and we find 

I min [X A ] = 2a p + 2 = 2 In \2 ^\P 2 Q 2 -J 2 \^j + 2. (40) 
Substituting further this result in (|20l ) we finally obtain 



JZ>8ii\/\P 2 Q 2 -J 2 \. (41) 

The extremal stationary near horizon geometry is in fact defined by equation (|32l ), by the 
same boundary conditions and by the same class of functions as those in the variational 
problem. Therefore, it is clear that the equality is saturated only for the extremal stationary 
near horizon geometry. This completes the proof. 

Remark. The case P 2 Q 2 = J 2 is formally outside of the class of functions we consider. 
In the language of stationary solutions, it corresponds to an extremal (naked) singularity with 
zero area. 

It is interesting to note that when PQ = 0, but P 2 + Q 2 ^ 0, the lower bound of the area 
depends only on the angular momentum but not on the nonzero charge in contrast with the 
Einstein-Maxwell gravity. 



4 Area-angular momentum-charge inequality for dilaton 
coupling parameter < y 2 < 3 

As we mentioned above, finding of sharp lower bound for the area J? in the case of arbitrary 
y seems to be hopeless. However, an important estimate for the area can be found for dilaton 
coupling parameter satisfying < y 2 < 3. The result is given by the following 

Theorem. Let H be a smooth, spacetime stably outermost axisymmetric marginally outer 
trapped surface in a spacetime satisfying the Einstein-Maw ell- dilaton-matter equations ([/]) 
with a dilaton coupling parameter J, satisfying < J 2 < 3. If the matter energy -momentum 
tensor satisfies the dominant energy condition and the dilaton potential is non-negative, then 
for every J in the given range, the area of "S satisfies the inequality 



A > %TlsJ \Q 2 P 2 - J 2 \, (42) 

where Q, P and J are the electric charge, the magnetic charge and the angular momentum 
associated with ( B. The equality is saturated only for the extremal stationary near horizon 
geometry of the y 2 = 3 Einstein-Maxwell-dilaton gravity with V (cp) = and T a t> = 0. 
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Proof. Let us first focus on the case < y 2 < 3 and consider the metric 

dL 2 = G AB dX A dX B (43) 

dX 2 + (d% + 2<S>d y i , -m'd<i>) 2 e~ 2 Y W 2 + e 2 ?W 2 f , , 

= — = — H h — dqr 

AX 2 X 3 Y 

and the associated functional 

= /.,{^ + 1 + ( 1 -^"-T^}^ < 44 > 

It is easy to see that I[X A ] > I[X A ] and therefore 

^>47c e («- 2 V 2 . (45) 
Redefining now the scalar field (f> = -^cp we find that the functional I[X A ] reduces to the 
functional I[X A ] for the critical coupling y 2 = 3. Hence we conclude that 



%%\J\Q 2 P 2 -J 2 \ (46) 

for every y with < y 2 < 3. 

The case y = needs a separate investigation. Fortunately, it can be easily reduced to the 
pure Einstein-Maxwell case. Indeed, it is not difficult to see that for y = we have 

I[X A ] >I EM [X% (47) 
where I EM [X A ] is the functional for the pure Einstein-Maxwell gravity. In Einstein- 
Maxwell gravity it was proven in that A > 8n\ J 2 + \{Q 2 +P 2 ) 2 which gives A > 



8n^/J 2 + i(Q 2 +P 2 ) 2 > Sn^/\Q 2 P 2 -J 2 \. 

5 Discussion 

In the present paper we derived area-angular momentum-charge inequalities for stable 
marginally outer trapped surfaces in the four dimensional Einstein-Maxwell-dilaton theory 
for values of the dilaton coupling parameter less than or equal to the critical value. The cou- 
pling of the dilaton to the Maxwell field leads in general to inequalities that can be rather 
different from that in the Einstein-Maxwell gravity. Some estimates for the sector y 2 > 3 
could be found if we impose the additional condition on the dilaton potential to be convex. 
We leave this study for the future. 

Given the current interest in the higher dimensional gravity it is interesting to extend the 
area-angular momentum-charge inequalities to higher dimensions. This is almost straight- 
forward in the case of Einstein equations ll22l . However, in the case of Einstein-Maxwell 
and Einstein-Maxwell-dilaton gravity the extensions of the inequalities is difficult. The cen- 
tral reason behind that is the fact that even the stationary axisymmetric Einstein-Maxwell 
equations are not integrable in higher dimensions [23J. Nevertheless, some progress can be 
made and the results will be presented elsewhere [[24]] . 

Acknowledgements: This work was partially supported by the Bulgarian National Sci- 
ence Fund under Grants DMU-03/6, and by Sofia University Research Fund under Grant 
148/2012. 



10 



References 

[1] A. Acena, S. Dain and M.E. Gabach Clement, Class. Quant. Grav. 28 105014 (2011); 
llarXiv: 1012.241 3r gr-qcll . 

[2] S. Dain and M. Reins. Phys. Rev. Lett. 107, 051 101 (201 1); HarXiv:1102.5215r gr-qcll. 

[3] M.E. Gabach Clement, HarXiv:1102.3834r gr-qcll. 

[4] J. L. Jaramillo, M. Reiris and S. Dain, Phys. Rev. D84, 121503 (2011); 
HarXiv:1106.3743r gr-qcll. 

[5] S. Dain, J. L. Jaramillo and M. Reiris, Class. Quantum Grav. 29, 035013 (2012); 
HarXiv:1109.5602r gr-qcll. 

[6] M.E. Gabach Clement and J.L. Jaramillo, llarXiv: 1111 .6248r gr-qcll . 

[7] W. Simon. Class. Quant. Grav. 29, 062001 (2012); llarXiv: 1109.6140r gr-qcl.l 

[8] M. Ansorg, J. Hennig and C. Cederbaum. Gen. Rel. Grav. 43, 1205 (2011); 
llarXiv: 1005.3 128l gr-qcll. 

[9] J. Hennig, M. Ansorg and C. Cederbaum. Class. Quantum Grav. 25 162002 
(2008): HarXiv:0805.43 20rgr-qcll. 

[10] J. Hennig, C. Cederbaum and M. Ansorg, Commun. Math. Phys. 293, 449 (2010); 
HarXiv:0812.281H gr-qcl ] . 

[11] S. Dain. Class. Quant. Grav. 29, 073001 (2012), HarXiv:1111.3615r gr-qcll. 

[12] G. Gibbons and K. Maeda, Nucl. Phys. B298, 741 (1988); 

[13] D. Garfinkle, G. Horowitz and A. Strominger, Phys. Rev. D43, 3140 (1991); D45, 
3888, 1992 (E). 

[14] D. Maison, Gen. Rel. and Grav. 10, 717 (1979). 

[15] S. Hojman, M. Rosenbaum and M. Ryan, Phys. Rev. D17, 3141 (1978). 

[16] S. Yazadjiev, Phys. Rev. D82, 124050 (2010); llarXiv: 1009.2442r hep-th11 

[17] B. Kleihaus, J. Kunz and F. Navarro-Lerida, Phys. Rev. D69, 081501 (2004); 
[arXiv:0309082[gr-qc]]. 

[18] S. Hayward, Phys. Rev. D49, 6467 (1994). 

[19] L. Andersson, M. Mars and W. Simon. Phys. Rev. Lett. 95, 111102 
(2005);[arXiv:0506013[gr-qc].] 

[20] L. Andersson, M. Mars and W. Simon. Adv. Theor. Math. Phys., 12(4), 853 (2008); 
HarXiv:0704.28 89rgr-qc11. 



11 



[21] I. Booth and S. Fairhurst, Phys. Rev. D77, 084005 (2008); [arXiv:0708.2209[gr-qc]]. 
[22] S. Hollands, Class. Quant. Grav. 29, 065006 (2012); larXiv: 11 10.58 M gr-qcll. 
[23] S. Yazadjiev, JHEP 1106, 083 (2011); llarXiv: 1 104.03781 Thep-thll. 
[24] S. Yazadjiev, In preparation 



12 



